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Abstract 

The results of a recent paper [0808.2909] are generalized. A more detailed 
proof is presented that under essentially all conditions, the non-linear clas- 
sical equations governing matter and gravitation in cosmology have "adi- 
abatic" solutions in which, far outside the horizon, in a suitable gauge, 
the reduced spatial metric gjj(x, t)/a^(t) becomes a time-independent func- 
tion Qij(x), and all perturbations to the other metric components and to 
all matter variables vanish. The corrections are of order a~^, and their x- 
dependence is now explicitly given in terms of Qij{x.) and its derivatives. The 
previous results for the time-dependence of the corrections to gij (x, t) / a'^ (t) 
in the case of multi-scalar field theories are now shown to apply for any 
theory whose anisotropic inertia vanishes to order a~^. Further, it is shown 
that the adiabatic solutions are attractive as a becomes large for the case of 
single field inflation and now also for thermal equilibrium with no non-zero 
conserved quantities, and the 0(a~^) corrections to the other dynamical 
variables are explicitly calculated in both cases. 
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I. Introduction 



In a recent paper [1], I addressed the problem of the evolution of non- 
Gaussian correlations during the long period when the observationally in- 
teresting wave lengths were outside the horizon, and the universe passed 
through transitions such as reheating and baryonsynthesis that are not well 
understood. It was shown on general grounds that, apart from quantum cor- 
rections, in a suitable coordinate system, the exact non- linear field equations 
always have an adiabatic solution for which Qij/a^ becomes time independent 
outside the horizon, and the perturbations to 300 and giQ and to all matter 
variables become negligible. This, however, is only part of the problem. In 
order to conclude in some class of theories that cosmological perturbations 
were really described by the adiabatic solution throughout the time when 
they were outside the horizon, it is necessary to show that for these the- 
ories the adiabatic solution is attractive in the limit of large a{t). In [1] 
the corrections to this solution for the metric were calculated explicitly for 
a theory with any number of scalar fields and an arbitrary potential, and 
it was shown that for all such theories this solution is attractive as far as 
the spatial part of the metric is concerned, but this could be shown for the 
scalar field perturbations only in the case of single field inflation [2]. The 
present paper clarifies these general arguments in Section II, and shows in 
Section III that the explicit solution for the metric obtained in [1] applies to 
any theory that (like scalar field theories) has vanishing anisotropic inertia 
to lowest order in perturbations. Sections IV- VI present an an improved 
argument that this solution is attractive for all metric and matter variables, 
not only in the case of single field inflation, but also for a sufficiently long 
period of local thermal equilibrium with no non-zero conserved quantities. 

II. The General Adiabatic Solution 

In this section we will give a general broken-symmetry argument that 
shows the existence of certain adiabatic solutions of the non-linear field 
equations for the metric and matter variables. The discussion will follow 
along the same lines as in Section II of [1], but we will here obtain more 
detailed results, that are of some interest in themselves, and will be used in 
the next section. 

Before proceeding to this task, it should be emphasized that the univer- 
sal existence of adiabatic solutions of the field equations does not in itself 
mean that these are the solutions that apply to the fluctuations in matter 
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and gravitation actually present in our universe. We will be able to show in 
Sections IV- VI that these solutions are attractive in two cases: single field 
inflation, and local thermal equilibrium with no non-zero conserved quanti- 
ties. Based on our experience with the linearized equations, it is plausible 
that fluctuations in these cases are in the basin of attraction of the adiabatic 
solutions, but we will not try to prove this for the non-linear equations. 

Turning to the existence of the adiabatic solutions, our discussion in this 
section is based on two very general assumptions: 

1. We assume that the dynamical equations governing matter and grav- 
itation are generally covariant. This of course is just a way of imple- 
menting the Principle of Equivalence. 

2. We assume that whatever the dynamical equations governing the met- 
ric and matter (including radiation) variables may be, these equations 
have a solution in which the metric takes the flat-space Robertson- 
Walker form, with goo = —1, goi = 0, and gij = gij/a? = 6ij, and 
in which all matter variables take their unperturbed form; that is, all 
densities and pressures and scalar fields are functions only of time, and 
all velocities and other 3-vectors vanish. 

We will show that under these conditions, for a suitable choice of space- 
time coordinates, these equations always also have a family of solutions that 
for large a{t) have the following properties: 

1. The metric for any of these solutions has components with 



2. Whether or not the energy and momentum of any particular con- 
stituent of the universe is separately conserved, for this solution its 
energy-momentum tensor has the form 



giji^, t) = a\t) [g^J{^) + h{t)nij{^) + /2(^)a^j(x)7^(x) + . . . 

5io(x,t) = /3(t)9,7e(x) + ... , 
5oo(x,i) = -1 + /4(^)7^(x) + . . . , 



(1) 



rio(x,t) = 53(t)Wx) + ... , 

roo(x,t) =p(t)+54(Wx) + ... , 



(2) 
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3. Four-scalars s{x,t) such as temperatures, number densities, or scalar 
fields, have the form 



s(x, t) = s{t) + h{t)n{x) + ... . (3) 

Here ^ij(x) is an arbitrary positive matrix function only of the spatial co- 
ordinates, which can conveniently be taken as the value of gij{x.,t)/a^{t) 
at some time t = T; TZij{x.) and TZ{x.) are the Ricci tensor and curvature 
scalar for the metric ^jj(x); the functions fn(t), Qnit) and h{t) are of or- 
der a~'^{t); dots denote terms of higher order in l/a{t); and a bar over any 
quantity indicates its unperturbed value. Although these solutions exist for 
any Qij (x) , quantum fluctuations in inflation produce a particular stochastic 
^jj(x), whose properties will not concern us here. Aside from quantum loop 
effects, a term of order a~" makes a contribution to correlation functions 
with characteristic wave number k that is expected to be suppressed outside 
the horizon by factors of order (k/aH)'"", where as usual H = a/ a. 

To illustrate what is meant by functions of order a~^{t), we note here 
that, as we will show in the next section, in suitable coordinate system, as 
long as there is no anisotropic inertia to order a~^, the functions fi{t) and 
/2(t) are in general given by 

/■* dt' /■*' 

/i(*) = 2/ ait")dt" (4) 



It a^t') Jt 

f2(t) = — / -^T—-- a(t")dt" + - -. — ^-^ / a(t")dt" 

2Jt a?{t') Jt ^ ' 2Jt H{t')a?{t') Jt ^ ' 

1 /•* H{t')dt' 



2Jt a'^{t')H{t') ' 



(5) 



while the coordinate system is defined so that /3(t) = 0. Here T is an 
arbitrary time, at which we choose gij/a^ to equal Qij. For instance, during 
the radiation-dominated era when a oc t^/"^ , at a time t^T, these give 

h{t) ^ . (7) 

The space derivatives in 7^jj(x) contribute two factors of a characteristic 
wave number k, so that these 0(a~^) terms do indeed make contributions 
of order (k/aH)'^. 
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On the other hand, the function f4{t) depends on the nature of the 
constituents of the universe. We show in Section V that during single field 
inflation 

while in Section VI we show that in a period of thermal equilibrium with no 
non-zero conserved quantities, 

-'^^ ' H^{t) Jt \H{t') 6H^{t') J V «^(* ) «^(* ) Jt ^ ' ) 

- (-^ + ^ t a{t') dt') . (9) 

2H{t) V a^{t) a3(i) Jt ^ ' ) ^ ' 

The first step in proving the existence of these solutions is to consider 
instead a trial configuration, in which Qijix, t) = a'^ {t)Qij (x) , where 0ij{x) is 
an arbitrary positive matrix function of position, and all other variables are 
the same as for the Robertson-Walker solution; that is, goo = —1, gio = 0; 
all densities and pressures take their unperturbed values; and all velocities 
vanish. Note that this is not a small perturbation to the Robertson- Walker 
solution, for no assumption will be made that Qij{x) is close to 6ij. If Gij{x) 
were constant then this would be an exact solution, since it could be ob- 
tained from the Robertson-Walker solution by a coordinate transformation 
X* — > A^jX^, where A^iA^j = Qij. Thus we expect the trial configuration 
to fail to satisfy the field equations only by terms involving derivatives of 
Qij{x.). These terms must be accompanied with factors of 1/a, as required 
by the condition that these field equations must be invariant under the scale 
transformations 

x'^\x\ a{t)^a{t)/X, (10) 

under which gij/a^, gio/a, goo, Tij/a^, Tio/a, Too (for each constituent of 
the universe), and all pressures, densities, temperature, etc. are invariant. 
(This is just a different way of expressing invariance under the coordinate 
transformations Ax*, under which g^y and T^i, transform as tensors, 

and a{t) is invariant.) Hence for sufficiently large a(i), it is a reasonable 
ansatz to seek a solution in which deviations of all scale-invariant metric 
and matter variables from the trial configuration are small perturbations. 

The field equations for the complete set of scale-invariant perturbations 
then take a form that can be symbolized as Dq = S, where (7 is a column of 
scale-invariant perturbations, D is a matrix formed from a linear combina- 
tion of first and higher time derivatives whose coefficients are scale-invariant 



5 



x-independent functions of time, and iS is a column of source terms, given 
by a linear combination of appropriate x-dependent 3-tensors formed from 
derivatives of ^jj(x), with factors of a{t) as required to make the source terms 
scale-invariant. Concrete examples will be given in Sections III, V, and VI. 
But even without going into the details of these equations, we can see the 
general features of a class of solutions. If we make the ansatz that the solu- 
tion for each perturbation is a linear combination of the same x-dependent 
3-tensors that appear in the source terms for this perturbation, with coeffi- 
cients given by functions only of time, then the dynamical equations reduce 
to a set of coupled ordinary inhomogeneous differential equations for these 
coefficient functions, with source terms having whatever powers of l/a{t) 
were required by scale invariance in the original equations. Such equations 
will always have an adiabatic solution, in which the functions of time ac- 
companying each 3-tensor are of an order in l/a(t) that is the same as the 
number of powers of l/a{t) in the corresponding source term, plus some 
number of sohitions of the homogeneous equations Dq = that may in gen- 
eral increase or decrease as a{t) increases. It is the time-dependence of the 
solutions of the homogeneous equation that tells us whether the adiabatic 
solution is attractive or not, the question addressed in Sections IV through 
VI. 

For gij/a? the source term must be a linear combination of the scale- 
invariant tensors o~^(t)7?.ij(x) and a~^(t)^ij(x)7?,(x), plus terms with more 
than two derivatives that are suppressed by more than two powers of l/a{t). 
For gio/a the source term must be a term proportional to the scale-invariant 
3-vector a~^{t)diTZ{x), plus terms with more than three derivatives that are 
suppressed by more than three powers of l/a{t). For goo the source term 
must be a term proportional to the scale- invariant scalar a~^(t)7?.(x), plus 
terms with more than two derivatives that are suppressed by more than 
two powers of l/a{t). Thus we expect these equations to have solutions 
of the form (1). The same reasoning applies to whatever matter variables 
(densities, pressures, velocities, and if necessary anisotropic inertia) that 
enter in the energy-momentum tensors for each constituent of the universe, 
giving energy- momentum tensors of the form (2) and scalars of the form (3). 

It should be noted that Eq. (2) tells us that STio is a gradient to order a~^, 
so that to this order there is no vorticity. On the other hand, Eq. (2) allows 
terms in STi j that arc not proportional to 6ij, so there is no general argument 
against the presence of anisotropic inertia to order a^^. The results of [1] 
show that there is no anisotropic inertia to this order in a theory of multiple 
scalar fields with an arbitrary potential, but solution of the coUisionless 
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Boltzmann equation [3] shows that anisotropic inertia of order can be 
produced in the Hnear approximation by gravitational perturbations of this 
order in a gas of non-interacting relativistic particles. 

III. Explicit Solution for the Metric 

Wc will now verify by detailed calculation that, under the assumptions 
of the previous section, for large a{t) there is indeed a solution of the gravi- 
tational field equations for which, in a suitable spacetime coordinate system, 
the metric takes the form (1), and wc will calculate the functions fnit) ap- 
pearing in this solution. For the metric, we use the ADM parameterization 
[4]: 

500 = -iV' + 9ijN'N^ , goi = gijN^ = 

gOO = _JV-2 , giO = Ni/N^ ^ gij ^ {3)gij _ ^i^jj^^ ^ (n) 

where ^^^g^i is the reciprocal of the 3 x 3-matrix g^j. It will be convenient 
also to write 

9ij{^,t) = a^{t)gij{y.,t) , (12) 

where a{t) is the Robertson- Walker scale factor appearing in the unper- 
turbed solution. In this notation, and in units with SttG = 1, the gravita- 
tional field equations are 

yi[&j-5'jC\)=-NT''j , (13) 
-Q-'^g'^Rij - C'jC^i + {C\f = 2iV2T°° , (14) 

+ N^VkCij + V.V.A^) = -Tij + ^a'-g^.T^x , (15) 

where 5*-'(x, t) is the reciprocal of the matrix gij{x, t); Rij{x, t) is the three- 
dimensional Ricci tensor for the metric gij (x, t) ; and C^j (x, t) is the extrinsic 
curvature of the surfaces of fixed time 

&j = a-^g"' Ckj , Cki = ^ [2aagki + a%, - V^iVi - V.iVfe] , (16) 

where Vj is the three-dimensional covariant derivative calculated with the 
three-metric g^j. 
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In accordance with the remarks of the previous section, we tentatively 
assume a solution in which g^^j and ^oo + 1 are small perturbations for large 
a{t), of order l/a^(f). We also adopt a definition of space coordinates for 
which iV* = 0, so that gio = 0. We then have 

&j = HS'j + ej, (17) 
where is, like and SN, a quantity whose leading term is of order a~^: 



Co = \t [fffe,- - 2il SN -gui + 0(a-^) . (18) 



Also in accordance with what we found in the previous section, the total 
energy momentum tensor for this solution takes the form 



Tij = a pcjij + 5Tij , Tio = 6TiQ Tqq = p + 6Tqo , (19) 

where p = 5H^, p = ~2H — 3H^, and STij/a'^, STiQ, and 6Toq are all of order 
a~^. The gravitational field equations now read 

V^ (6 - S'j^'^k) = 5Toj + 0(a-^) , (20) 

-12H^SN + 25Tqq = -a-^g'^Rij + AH^\ + 0(0"^) , (21) 

Cj + 3HCj + H6'j^''k + {^H^ - H)SNS'j = a-'^g'''Rkj 

+ a-^g^'STkj - ^a-^S^jg^'STki + ^S^jSToo + 0{a-^) . (22) 

Now, as remarked in Section II, to order a^^ these solutions have no 
vorticity, so we can introduce a momentum potential U, of order a~^, such 
that 

6Tio = diU . (23) 
The equation of momentum conservation then reads 

a-^g'^Vi5Tkj = 2Hdj5N + dj{U + mU) + 0(0"^) . (24) 

We can always also write 

5Tij=a'^(gij5p + Xlij) , (25) 

where 11^^ is a 3-tensor of order a^^ representing anisotropic inertia, with 
g'^^Iiij = 0. Eq. (2) for the total energy momentum tensor together with 



8 



the Bianchi identity for TZij tells us that g'^^VidT^j is a gradient, so we can 
introduce another potential 11 such that 

g'''ViUkj = djU . (26) 

The momentum-conservation equation (24) then reads 

n = -Sp + 2HSN + U + 3HU . (27) 

Using Eqs. (21), (25), and (27), the field equation (22) now takes the form 



1 



1 



r.kl ; 



where 



+ ~g'^^k3 + \5'o^ + 0{a-^), 



t« . I .TJ 
" J — S J 2 



(28) 



(29) 



As a check, we note that the field equation (20) may be put in the form 

Vi - (5^,S\) + 0(a-4) = , (30) 

which is automatically consistent with Eq. (28) as a consequence of Eq. (26) 
and the Bianchi identity obeyed by Rij. 

Up to this point we have kept open the possibility of anisotropic inertia 
terms of order in order to emphasize the very great simplicities offered 
if anisotropic inertia can be neglected to this order. If we now assume that 
Iljj = to order a~^, so that also 11 = to this order, then Eq. (28) reads 



S'-; + 3i/S*o — —77 



+ 0{a-^) . 



(31) 



The general solution of Eq. (31) is of the form 



2Mx,t) 



6?^'=(x)7^fe,•(x) - -5^,•g^'(x)7^fe^(x) 



1 



a^{t) Jt 



a{t') dt' 



(32) 



where T is any fixed time, ^ij(x) is the value of gij{x, t) at that time, 7^ij(x) 
is the Ricci tensor calculated from the 3-metric ^ij(x); and -B'j(x) is some 
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function of x (and T), appearing in the solution of the homogeneous equation 
corresponding to Eq. (31). It is striking that once we assume an absence of 
vorticity and anisotropic inertia, all other terms that depend on the details 
of the energy-momentum tensor cancel in this solution. 

To solve for the metric, we need to complete our choice of gauge. By 
using Eqs. (18), (21), and (29), we have 

IH^ H 

tij = 2gikE''j + -^gijE\ - —-^hfg'^^Rkl + hjX + 0{a-^) , (33) 
where X = 0(l/a^) depends on the matter perturbations 

X = -| (,5roo - 2(3ij2 + H)5n) -ui^ + ^ j . (34) 

(There arc other ways of writing Eq. (33), with apparently simpler definitions 
of X, but as we will see the definition used here provides special advantages.) 
Under a shift t ^ t + e(x, t) in the time coordinate, with e of order 1/a^ 
(and a corresponding transformation — > + Q'^^ J dta'''^de/dx^ to keep 
Ni = 0), the quantity X undergoes the transformation 

X^X + 2-{eH)+... , (35) 

so we can evidently choose e to make X = 0. This choice is not unique; 
after choosing space and time coordinates so that giQ = and X = 0, we 
can preserve both conditions by making a further gauge transformation with 

.^t + r(x)/fi(,), + (36) 

where r is an arbitrary function only of x. This residual gauge freedom will 
be important in Sections V and VI. 

In a theory of multiple scalar fields X is a linear combination of terms 
proportional to the scalar field perturbations and their time derivatives, so 
the choice X = is a generalization of the choice of gauge commonly used 
in studying single field inflation, in which the scalar field is unperturbed. In 
general, the gauge choice X = provides the great advantage that we can 
solve Eq. (33) for gij without worrying about the matter variables to which 
the metric is coupled: 



gij{^,t) = ^ij(x) + 2 



1 1 /■* fit' /■*' 

T^iji^) - -a.,(x)g^'(x)7^feKx)J ^ a{t")dt" 
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ft Jf' 
+ 2^.(x)i.^(x)/ 




H{t')dt' 

't a\t')H{t') 



.11 




H^{t')dt' 



+ 0(a-\t)) , 



(37) 



where T is again any fixed time, and Qiji'x.) and 7?.ij(x) are the values of 
gij and the associated Ricci tensor at that time. This is the same solution 
found in [1] for the special case of multiple scalar fields with an arbitrary 
potential. In the language of Eq. (1), this solution tells us that 



Of course, having chosen space coordinates so that gio = 0, we have /^{t) = 
0. The calculation of f4{t) is taken up in Sections V and VI. 



We have seen in Eq. (34) that the general solution of the dynamical 
equations for matter and metric variables that differs from the "trial con- 
figuration" (for which g^j {x, t) / (t) = ^jj(x), giQ = 0, 500 = ~1) and all 
matter variables are unperturbed) by small terms, of order a""^, has a spa- 
tial metric that approaches Gij{x) for large a{t). But in order to show that 
the trial configuration is truly attractive, we also need to consider the other 
metric component 500 = — -/V^ (the coordinate system having been chosen 
so that gio = 0), and the various matter variables. We will do this in the 
next two sections for two cosmological models: single field inflation and local 
thermal equilibrium with no non-zero conserved quantities. Both of these 




(38) 



(39) 



IV. Other Variables 
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models have vanishing anisotropic inertia, and therefore have a spatial met- 
ric described by the results of the previous section. Prom the perfect-fluid 
form of the energy-momentum tensor, we have in both cases 

5Too = 5p + 6H^6N + 0(a"^) (40) 

so Eqs. (21) and (29) give 

26p = -a'^g'^Rij + AHE'^k - GHU + 0(a"'') , (41) 

while Eqs. (40) and (34) show that the gauge condition X = gives 



= - 

H 



(Sp - 2H6N^ -u(^l + (42) 



Finally, for zero anisotropic inertia the equation (27) of momentum conser- 
vation gives 

= -dp + 2HSN + 3HU + U + 0{a-^) . (43) 

This gives three relations among the four quantities 5N, 6p, 6p, and U, so 
with one additional relation we can calculate all four. 

V. Single Field Inflation 

In general, for a single scalar field (p{t) + 5(/?(x, t), with a conventional 
kinematic term, the quantity X defined by Eq. (34) is given by [1]: 

H 



X = —(pS(p + -j^ (^(pS(p — 'pS(p'^ + 0{a ^) 



1 



-2HS(p - ^ ( - H6(p + 2Hd(pj 



V-2H 

The gauge choice that X = to order tells us then that 



+ 0{a-^) . (44) 



Mx,i) = /(x)^^-^^, (45) 

where /(x) is some function only of x. Under a residual gauge transforma- 
tion of the form (36) (which preserves the conditions X = and gio = 0) 
the scalar field perturbation is shifted by 

A(5(^(x,0 = -^{t)T{x)/H{t) = -^-2H{t)r{^)/H{t) , 
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so by choosing t(x) = /(x), we can make 5(p = 0, completing our choice of 
gauge. The momentum potential defined by Eq. (23) is here 

U = 06ip , (46) 

so the gauge choice S(p = also entails U = 0. The density and pressure 
perturbations here are 

5p = 5p = ^6goo^^ = 2HSN , (47) 

so Eqs. (42) and (43) are satisfied, and the only variable that needs to be 
examined to check the attractive nature of the adiabatic solution is c^^oo = 
-26N. It is given by Eqs. (47) and (41) as 



SN = ^ 
AH 



AH 



^ + ^ / a{t') dt' 



+ Oia-^) . (48) 



This is 0(o ^), and hence vanishes for large o, so the adiabatic solution is 
indeed an attractor in the case of single field inflation. In the language of 
Eq. (1), Eq. (48) gives 

In the gauge adopted here there are no matter perturbations during single 
field inflation, and so nothing else to check. 

VI. Local Thermal Equilibrium 

In thermal equilibrium with no non-zero conserved quantum numbers 
the pressure and energy density in any gauge are functions only of the tem- 
perature, so 

dp = {p/p) bp = ^-1 - ^ j Sp . (50) 

Using this and Eqs. (41) and (42) in Eq. (43) gives a differential equation 
for the momentum potential U : 



13 



where T may be taken as any time during the period of thermal equihbrium, 
most conveniently at its beginning, and Qij and TZij are the reduced metric 
Qij/a? and associated Ricci tensor at that time. The solution is 

C(x,t) = /(x)f| 



X 



H{t)JT \H{t'))\ 6H{t')H{t')J 



where /(x) is an arbitrary function of position. At first sight, this does 
not appear very attractive, in either the mathematical or the colloquial 
sense. The first term, representing a solution of the homogeneous equation 
corresponding to (51), is of zeroth order in 1/a, and so does not become 
small for large a. But this term is a gauge artifact. From the definition (23) 
of the momentum potential, we can easily see that under the residual gauge 
transformation (36) (which preserves the conditions X = and gio = 0), U 
transforms as 

C/(x,t)^C/(x,t)+(^^^r(x) (53) 

so by choosing the arbitrary function r(x) to have the value — /(x)/2, we 
can cancel the first term in Eq. (52). The remainder is of order a~^, and 
hence vanishes for large a. 

The remaining perturbations 6goo = —25N and 6p are algebraically re- 
lated to U by Eqs. (41) and (42), which give 

(5p(x,t) = -SH{t)U{^,t) + lg^^■(x)7^,,■(x) (--i^ + ^ £dt' a{t': 

(54) 
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and 

2mm^.t) = (11) t/(x.t)+i5«(x)K,(x) 

(55) 

Both are of order a ^, so in this case, too, the adiabatic solution is com- 
pletely attractive for large a. In the language of Eq. (1), Eqs. (55) and (52) 
give: 

■'^^^ H^{t) Jt \H{t') QH^{t')j\ a^{t') a^{t') Jt ^ ' J 
1 / 1 H{t) 



2H{t) V a2(t) a3(i) 



+ fix (^"^ 
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